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INVESTIGATION  OF  EXPLOSION  GENERATED  SV  Eg  WAVES  IN  2-D 
HETEROGENEOUS  CRUSTAL  MODELS  BY  FINITE-DIFFERENCE  METHOD 

SUMMARY 

A  linear  finite-difterence  (FD)  method  was  used  to  compare  the  excitation  of  far-held  !’- 
and  SV-waves  generated  by  shallow  dilatational  sources  in  a  suite  of  heterogeneous  2-D  crustal 
models.  The  crustal  models  tested  included  simple  layered  structures,  media  with  random  velo¬ 
city  perturbations  having  Gaussian  or  self-similar  autocorrelation  functions,  media  with  rough 
or  gentle  topography  generated  by  Markov  chains,  and  laminated  media  with  sinusoidal  folds. 
The  numerical  experiments  were  conducted  by  directing  a  broadband  planar  P-  or  SV- wave 
with  appropriate  incidence  angle  upon  the  testing  models.  The  dilatational  strain  history  at  a 
shallow  linear  array  of  grid  points  was  then  recorded  so  that  the  far-held  P-  or  SV(Lg)-waves 
from  shallow  dilatational  sources  could  be  inferred  by  use  of  the  principle  of  reciprocity.  The 
raw  FD  synthetics  were  deconvolved  so  as  to  represent  the  response  due  to  explosion  sources 
with  a  fixed  yield.  The  mean  peak  amplitude  of  the  synthetics  for  each  model  are  compared  to 
that  for  a  reference  model  consisting  of  a  simple  layered  medium.  The  average  energy  content 
in  an  appropriate  signal  window  was  measured  as  a  complement  to  the  amplitude  measurement. 
Both  approaches  show  essentially  the  same  pattern  of  P/SV  excitation,  namely  that  models  with 
topography  consistently  produce  the  strongest  P-SV  conversion  among  all  types  of  crustal 
models.  The  introduction  of  interfaces  (e.g.,  dipping  layers)  alone  does  not  by  itself  increase 
SV  excitation  with  the  required  slowness  range.  Thus  mh(P)  •  mh(f  g)  appears  to  be  smaller  for 
models  with  topographic  relief  (e.g  ,  the  Novaya  Zcmlya  Island)  than  for  models  with  dipping 
layers  or  folded  sedimentary  rocks  (e.g  ,  Shagan  River,  eastern  FKTS). 

These  synthetic  results  are  consistent  with  observations  for  Novaya  Zemlya  (Nuttli,  1988) 
and  Shagan  River  (Nuttli, 1986),  based  on  WWSSN  film  chip  readings  of  Lg.  The  Novaya 
Zemlya,  which  has  rough  topography,  has  a  higher  Lg  relative  to  P  (  mh(P)  -  mh(Lg)  =  0.04  ) 
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than  does  the  somewhat  flatter  Shagan  River  test  site  (  mb(P)  -  mb(Lg )  =  0.27  ).  However, 
Nuttli  (1987)  also  obtained  relatively  low  mb(Lg)  for  the  Degelen  Test  Site,  which  is  only  70 
km  away  from  Shagan  River.  If  this  Degelen-Shagan  River  bias  ia  real,  it  is  not  explained  by 
the  FD  results  obtained  to  date.  However,  some  recently  archived  high-quality  digital  seismo¬ 
grams  recorded  at  the  Chinese  Digital  Seismic  Network  indicate  more  Lg  excitation  (with 
respect  to  P)  at  Degelen  than  at  Shagan  River,  which  is  consistent  with  the  numerical  rest  its. 
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INTRODUCTION 

The  seismic  Lg  wave  is  one  of  a  number  of  regional  phases  that  propagate  in  the  con¬ 
tinental  lithosphere.  Because  the  anclastic  attenuation  of  1-s  period  Lg  wave  is  small  in  shield 
and  geologically  old  stable  regions,  Lg  wave  amplitudes  provide  a  useful  tool  for  estimating 
1-s  period  magnitude,  such  as  mh  for  small  earthquakes  and  explosions  (Nuttli,  1973).  Further¬ 
more,  the  radiation  of  Lg  is  more  isotropic  than  that  of  P  and  S  waves,  which  adds  to  its  use¬ 
fulness  as  a  magnitude  estimator  for  small  events,  because  full  azimuthal  coverage  is  not  essen¬ 
tial  and  thus  reliable  magnitude  determin  >tiops  can  be  made  from  the  data  of  only  a  few  sta¬ 
tions  (Nuttli,  1986a).  Following  this  line,  Nuttli  (1986a, b,  1987,  and  1 98S)  compared  mh(P)  - 
mh( Lg)  at  the  Novaya  Zemlya,  Nevada  Test  Site,  the  Shagan  River  (eastern  portion  of  Fast 
Kazakhstan  Test  Site)  as  well  as  the  Degelen  region  (central  portion  of  Hast  Kazakhstan  Test 
Site)  as  tallows: 


TABLE  1.  mb( P)  - 

mh{ Lg)  AT  VARIOUS  SITES 

Site 

"<*,( P)  -  Wfr(Lg) 

Reference 

Novaya  Zemlya 

-0. 1 1±0.02 

Nuttli  (1988) 

Shagen  River  EEKTS 

-t-0.04t0.12 

Nuttli  (1986b) 

Degelen  CEKTS 

0.2/±0.03 

*T  •  '  r.?mr's 

i>uiui  ^  :  J\j  i  j 

N;  vada  Test  Site 

-0.31  ±0.02 

Nuttli  (1986a) 

The  large  difference  between  mh( P)  -  mb( Lg)  values  at  Shagan  River  and  Degelen,  0.23 
magnitude  unit,  is  quite  surprising  inasmuch  as  the  two  test  sites  are  less  than  70  km  apart. 
Without  knowing  the  actual  values  of  the  explosion  yields  at  the  Eastern  Kazakhstan  sites,  it  is 
not  obvious  at  all  to  decide  whether  the  mb( P)  or  the  mfc(Lg)  values  at  Degelen  are  anomalous. 
The  difference  of  0.23  magnitude  unit  implies  that,  for  an  explosion  of  given  yield,  either  the 
P-wave  amplitudes  at  Degelen  are  1.9  times  iarger  than  at  Shagan  or  that  the  Lg-wave  ampli¬ 
tudes  are  1.9  times  smaller  at  Degelen.  Of  course,  both  the  P  and  Lg  mioht  be  different 
between  Degelen  and  Shagan  for  a  given  yield. 
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A  number  of  possible  explanations  of  the  differences  between  mh( P)  -  m^L.g)  values  at 
Shagan  River  and  Degelen  have  been  discussed  (Nuttli,  personal  communication).  The  first  is  a 
difference  in  coupling  of  P  and  Lg  waves  at  the  two  sites.  Either  the  P-'vave  coupling  is 
significantly  more  efficient  at  Degelen  than  at  Snagan,  or  the  Lg-wave  coupling  at  Degelen  is 
less  than  at  Shagan.  A  second  possible  explanation,  somewhat  related  to  the  first,  is  that  mb( P) 
is  larger  at  Degelen  because  of  testing  practices.  Station  sampling  bias  might  be  considered  as 
the  third  explanation,  but  the  data  Nuttli  used  indicate  that  there  is  little  such  bias,  if  any.  A 
fourth  explanation  is  related  to  the  fact  that  most  of  the  Degelen  explosions  are  smaller  than 
those  at  Shagan,  and  thus  the  reported  mb{ P)  values  for  the  Degelen  events  might  be  overes¬ 
timated  because  only  stations  recording  large  amplitude  P  waves  were  used  to  obtain  the  Inter¬ 
national  Seismological  Center  mb(P)  values.  However,  Nuttli  (1987)  reported  that  the  largest 
Degelen  events  of  20  February  1975.  26  March  1978,  28  July  1978,  and  29  November  1978 
have  an  average  mh(P)  -  mb(Lg)  value  of  0.27,  the  same  value  as  foi  die  entire  set  of  Degelen 
exp’  sions.  The  fifth  possible  explanation  is  that  the  slope  of  mb(Lg)  versus  mb( P)  is  less  than 
unity  because  of  differences  in  comer  frequencies  in  spectra. 

It  is  the  purpose  of  this  report  to  check  the  first  hypothesis.  CEKTS  has  considerable 
topographic  relief  (Rodean,  1979)  whereas  EEK.TS  consists  of  folded  sedimentary  rnrks  (Nor- 
dyke,  1973).  We  therefore  investigated  several  2-D  models  for  the  excitation  of  SV  Lg  as  well 
as  P  waves  using  linear  elastic  finite-difference  calculations.  These  simplisuc  models  simulate 
die  excitation  of  far  field  SV  and  P  waves  from  dilatational  line  sources  embedded  in  2-D 
heterogeneous  models  with  and  without  free  surface  topography. 

Strictly  speaking,  no  adequate  analytical  theory  is  available  to  model  tire  near-source  cru¬ 
stal  effect  of  the  type  investigated  here.  Lg  is  thought  to  be  superpositions  of  higher  Rayleigh 
and  Love  modes.  In  terms  of  higher  modes,  the  complexities  of  mode  conversions  and  scatter- 
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2.  self-similar  media  generated  by  modulating  the  wavenumber-wavenumber  spectra  of  a 
white  velocity  held  with  the  2  D  Fourier  transform,  - - ,  of  a  special  Von  Karmon 

1  4 


correlation  function. 


3  iolded  lasers  of  sinusoidal  shape  with  specified  wavelength,  amplitude  and  velocity 

profile 

4.  dipping  sedimentary  layers  with  specified  velocity  profile. 

5.  simple  layered  models 

On  some  of  the  models  there  was  superimposed  fairly  tough  topography  generated  by  a 
Markov  chain.  Fouiteen  2-1)  geologic  structures  were  used  to  cab  1 1 1  ate  the  excitation  of  far- 
lield  SV  and  P  waves  from  ddational  line  sources.  Some  typical  simplistic  models  are  shown 
in  Figures  1  and  2.  Similar  random  media  have  been  used  widely  in  modeling  the  scattering  or 
propagation  of  either  acoustic  or  elastic  waves  by  the  liuite-diffetence  method  (Frankel  and 
Clayton,  1984,  1986;  l.evandei  and  Hill.  1985;  l.ecandcr.  1 1  > S 5 ;  l-'iankel  and  Wennerberg, 
1987;  Mcl  aughlin  and  Jib,  1987).  In  the  heterogeneous  portion  of  the  grid,  the  S-wave  velo¬ 
city  (P)  was  assumed  to  be  directly  related  to  the  P-wave  velocity  (it.)  by  the  following  linear 
relationship:  for  a  <  2.8,  (5  -  0.45  a;  for  3.2  <  a  <  4.5,  jf  .  0.50  <x;  for  a  >  5.2,  p  ■=  0.59  a. 
and  linear  interpolation  was  used  for  the  transition  intervals.  These  heterogeneous  media  are 
then  embedded  into  a  homogeneous  half  space  with  P-  and  S-wave  velocities  of  6  and  3.55 
kin; sec  respect i  vcl  v. 

I  he  incident  wave  is  a  broadband  Oluuka  planar  SV  wave  with  apparent  velocities  of  4  5 
km/sec  incident  upon  these  models  The  diluiion.il  strain  was  recorded  at  an  array  of  locations 
with  an  average  depth  of  0  5  km  in  the  grid.  By  use  t.f  the  reciprocity  principle,  the  displace¬ 
ment  response  at  lat  distance  (lelescismi  )  was  determined  t-u  a  dil.itiona!  Ime  source  with  a 
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von  Seggern  and  Blandtoid  (1972)  reduced  displacement  potential  lot  50KT  in  hard  rock.  The 
same  technique  has  been  utilized  previously  in  modeling  the  effects  of  die  geologic  structure  of 
Yucca  Valley  (McLaughlin  1 1  al ,  1987)  as  well  as  that  of  the  topographic  configura'ion  of 
Ahaggar  mountain  on  teleseismic  (McLaughlin  and  hh.  1986).  One  minor  difference  is 

that  here  the  incidence  unities  of  the  incoming  SV-  and  i’  wave  were  52"  and  20"  respectively, 
corresponding  to  the  appropriate  apparent  velocity  range.  To  assure  that  the  initial  incoming 
wave  lies  completely  in  the  homogeneous  ha  if  space,  we  have  selected  grid  dimensions  of  250 
by  ‘150  and  250  by  250  tor  S',  at?  **  wave  simulations,  respectively.  The  mesh  spacings 
were  Ax  -  Ar  0.050km.  Thus  the  two  gods  ate  of  the  same  width  12.5km  and  two  different 
depths:  21. 5km  and  17  5km.  l  ine  to  the  i  mu.ci  m  of  die  end.  a/e.  the  heterogeneous  crustal 
layers  were  set  to  2  oi  5  km.  This  dicmcm  a  with  the  realistic  crustal  depth  should  be  kept  in 
mind  when  interpreting  the  results  ■Mila  ugh  the  temporal  spacing  was  0.005  sec,  the  band 
beyond  5  Hz  is  less  reliable  due  to  the  Hie-.ite'de  grid  dispersion  of  finite-difference  simulation. 

The  FI)  Synthetics  were  deconvolved.' convolved  so  as  to  represent  the  outgoing  waves 
due  to  the  same  evpios.on  line  .sour. e  in  each  model.  The  results  from  these  simulati  in-  are 
pi ese tiled  in  Table  2  l  he  ex.  nation  oi  I’  .md  SV  waves  lot  each  model  are  referenced  to  the 
excitation  for  a  uniform  2  ktn  layer  (u  5.0km  s,  ?\  2.7-tkm/s)  over  a  uniform  half  space  (a 

«  O.Okm.s,  (5  5.55km.sj  Th  •  models  ate  art attged  in  order  of  decteasing  SV  Lg  excitation 

relative  to  the  teference  model  (j  0).  !)aia  shown  under  columns  P  and  Lg  are  log(7y/’o)  and 
/.cf  te-pcctivch .  where  both  L  and  |  g  Jie  t!ie  averaged  peak  amplitude  measurement 
o'  > ' w*  xvmlvtics  de:  tv  ed  I > t i '  .i  •convolute  m.  low  pass  lilteiane  <•■ 
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TABLE  2.  COMPARISON  OF  EXPLOSION  P  AND  SV  Lg  EXCITATION 


model 

P 

Lg 

P-Lg 

Description  of  the  model 

0 

0.000 

0.000 

0.000 

1 -uniform  layer  (5+0%,  2km  thick) 

1 

-0.207 

0.202 

-0.409 

rough  TOPO  +  1  uniform  layer  (5+0%,  2km  thick) 

2 

0.006 

0.132 

-0.139 

gentle  TOPO  +  self-similar  layer  (5+10%,  2km) 

3 

-0.196 

0.110 

-0.3C5 

rough  TOPO  +  Gaussian  layer  (5+10%,  2km) 

4 

-0.023 

0.073 

-0.096 

gentle  TOPO  +  1  uniform  layer  (5+0%,  2km) 

5 

-0.034 

0.044 

-0.078 

self-similar  layer  (5+10%,  2km  thick) 

6 

-0.162 

0.019 

-0.181 

folded  sinusoidal  layers(L=2, 11=2.5,5+20%) 

7 

-0.031 

0.014 

-0.045 

folded  sinusoidal  !ayers(L=2,H=2.5,5+I0%) 

8 

-0.134 

-0.037 

-0.098 

self-similar  layer  (5+20%,  2km  thick) 

9 

-0.029 

-0.037 

0.008 

folded  sinusoidal  layers(L=5,H=2.5,5+10%) 

10 

0.003 

-0.058 

0.061 

2-Gaussian  layer  (4.5+10%/5+ 10%,  total  2km) 

1 1 

0.019 

-0.091 

0.1 10 

steeply  dipping  layers  (52°) 

12 

0.011 

-0.093 

0.104 

gently  dipping  layers  (26n) 

13 

0.018 

-0.137 

0.155 

steeply  dipping  layers  (-52°) 

14 

0.009 

-0.143 

0.152 

gently  dipping  layers  (—26°) 

Some  observations  are  immediate: 

1.  For  self-similar  models,  weak  variation  in  the  medium  velocities  causes  more  Lg  genera¬ 
tion. 

2.  Dipping  layers  (models  1 1  through  14)  generate  smaller  Lg  than  the  normalizing  model. 

3.  Media  with  topography  (e  g.  models  1  through  4)  which  represent  CEKTS  all  generate 
more  Lg  than  the  normalizing  model. 

4.  Dipping  layers  (models  11  through  14)  are  more  efficient  than  all  other  models  for  P 
excitation. 

s  mh(P)  -  fti// Lg)  at  EEKTS  (e  g.  models  11  through  14)  are  larger  than  for  any  other 
model. 

We  have  noticed  that  the  smoothing  of  the  interlace  between  the  sediment  layer  and  the 
half  space  would  not  significantly  affect  the  result. 
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Observations  2,3,4  and  5  can  be  explained  easily  as  that  topography  creates  more 
apparent  explosion  P-SV(Lg)  coupling.  The  introduction  of  interfaces  alone  does  not  of  itself 
increase  SV  excitation  with  the  required  slowness  range.  Although  this  result  is  in  seeming 
contradiction  to  Nuttli’s  (1987)  CEKTS-EEKTS  observation,  our  previous  numerical  experi¬ 
ments  tend  to  support  this  explanation. 

As  a  comparison,  the  results  corresponding  to  frequency-domain  measurements  are  also 
presented  for  the  0. 5-1.0  Hz  bandwidth  in  Table  3.  Again,  the  average  spectral  level  of  the 
dilatational  strain  history  in  the  0. 5-1.0  11/  bandwidth  from  the  array  of  locations  is  determined 
relative  to  a  reference  model.  Note  that  the  results  under  the  Lg  and  P-Lg  columns  in  Table  3 
follow  the  same  pattern  as  the  time-domain  approach  with  the  exception  of  the  folded  models 
and  the  20%  self-similar  model. 


TABLE  3.  EXPLOSION  EXCITED  P  AND  SV  Lg  ON  0.5- 1.0  Hz 


model 

P 

LS 

P-Lg 

Description  of  the  model 

0 

0.000 

0.000 

reference  model 

1 

-0.400 

0.083 

-0.483 

rough  TOPO-t uniform  layer(5+0%,2km) 

2 

-0.049 

0.057 

-0.106 

gentle  TOPO+self-similar  layer(5  +  10%,2km) 

3 

-0.363 

0.063 

-0.426 

rough  TOPO+Gaussian  layer(5. 0+10%, 2km) 

4 

-0.180 

0.019 

-0.199 

gentle  TOPO+uniforni  layer(5+0%,2km) 

5 

0.016 

0.009 

0.007 

self-similar  layer(5+10%,2km) 

6 

0.099 

-0.031 

0.130 

folded  sinusoidal  layers(5+20%,L=2,H=2.5) 

7 

0.058 

-0.101 

0.159 

folded  sinusoidal  !ayers(5+10%,L=2,H=2.5) 

8 

-0.026 

-0.049 

0.023 

self-similar  !ayer(5+20%,2km) 

9 

0.015 

-0.163 

0.178 

folded  sinusoidal  layers(5+10%,L=5,H=2.5) 

10 

0.083 

-0.007 

0.090 

2-Gaussian  layer(4.5  +  10%/5. 0+10%, 2km) 

11 

-0.008 

-0.048 

0.040 

steeply  dipping  layers(52°) 

12 

-0.024 

-0.057 

0.033 

gently  dipping  layers(26”) 

13 

0.015 

-0.086 

0.101 

steeply  dipping  layers)— 52°) 

14 

-0.001 

-0.103 

0.102 

gently  dipping  layers(-26") 
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DISCUSSIONS  AND  CONCLUSIONS 

The  method  of  finite-difference  was  used  to  simulate  the  excitation  of  explosion  generated 
far-field  P  and  SV  Lg  waves  in  various  models  of  crustal  heterogeneity.  While  we  continue  to 
experiment  with  various  models,  our  preliminary  results  indicate  that  P  to  SV  conversion  is 
strongly  enhanced  by  velocity  variation  in  the  vicinity  of  rough  topography  and  the  introduc¬ 
tion  of  low  velocity  layers  near  the  surface.  The  introduction  of  interfaces  alone  does  not  of 
itself  increase  SV  excitation  with  the  required  slowness  range. 

These  synthetic  results  are  consistent  with  observations  for  Novaya  Zemlya  (Nuttli,  1988) 
and  Shagan  River  (Nuttli,  1986),  based  on  WWSSN  film  chip  readings  of  Lg.  The  Novaya 
Zemlya,  which  has  rough  topography,  has  a  higher  Lg  relative  to  P  (  mb(P)  -  mb(Lg)  =  0.04  ) 
than  does  the  somewhat  flatter  Shagan  River  test  site  (  mb(P)  -  >nb(Lg)  =  0.27  ).  However, 
Nuttli  (1987)  also  obtained  relatively  low  mb(Lg)  for  the  Degelen  Test  Site,  which  is  only  70 
km  away  from  Shagan  River.  If  this  Degelen-Shagan  River  bias  is  real,  it  is  not  explained  by 
the  FD  results  obtained  to  date.  However,  some  recently  archived  high-quality  digital  seismo¬ 
grams  recorded  at  the  Chinese  Digital  Seismic  Network  indicate  more  Lg  excitation  (with 
respect  to  P)  at  Degelen  than  at  Shagan  River,  which  is  consistent  with  our  numerical  results 
(Figure  41). 

Although  we  cannot  presently  explain  Nuttli’s  (1987)  Degelen  results,  we  predict  substan¬ 
tial  variations  in  SV  Lg  excitation  by  explosions  embedded  in  crustal  heterogeneity.  It  seems 
that  P-to-SV  is  not  the  only  mechanism  for  explosion  Lg  excitation,  so  it  is  necessary  to  inves¬ 
tigate  the  excitation  of  SH(Lg)  as  well. 

A  possibility  is  that  Nuttli’s  mb\ Lg  relationship  might  be  relaied  to  Rayleigh-to-P  conver¬ 
sion  away  from  the  immediate  location  of  the  source.  Our  numerical  simulations  treated  only 
ihe  P  S  conversions  that  might  occur  within  a  few  km  of  the  source,  and  given  some  simple 
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models.  If  either  the  Rayleigh  excitation  or  the  Rayleigh  scattering  is  different  for  CEKTS  and 
EEKTS,  then  we  could  see  the  difference  in  Rayleigh  to  Lg.  Since  the  two  locations  are  only 
70  km  apart,  the  Rayleigh-to-Lg  difference  would  presumably  have  to  occur  in  the  first  20-25 
seconds.  Thus  it  seems  necessary  to  examine  if  the  P-roda  are  different  for  Degelen  and 
Shagan  in  the  first  20-25  seconds.  Similarly,  P-SV  conversion  could  be  happening  further 
away  from  the  source  than  we  are  modeling,  ft  is  also  possible  that  the  non-linear  source 
effects  might  produce  larger  SV  at  one  site  versus  another.  These  hypothesis  as  well  as  the  3- 
dimensional  effects  were  not  addressed  in  our  current  experiments. 

Baumgardt  (personal  communication,  1987)  pointed  out  that  the  observations  at  NOR- 
SAR  for  Degelen  and  Shagan  events  (Ringdal,  1982)  did  not  show  the  bias,  and  he  proposed 
another  explanation  for  Nuttli’s  observations  as  that  Nuttli’s  Degelen  observations  were  made 
almost  entirely  with  stations  to  the  south,  whereas  the  Shagan  data  contain  observations  from 
the  NW  at  Scandinavian  stations.  His  station  corrections  for  Q  may  therefore  have  been  biased. 
Baumgardt  (1985)  proposed  that  Lg  loses  energy  as  it  crosses  the  Urals  on  its  way  to  the  Scan¬ 
dinavian  stations,  and  this  loss  may  not  be  accounted  for  in  Nuttli’s  Q.  It  is  interesting  to  note 
that  in  Nuttli’s  (1986b)  original  paper  on  Shagan  events,  his  revised  Q0  values  are  the  same  as 
or  less  than  the  original  Q0  values  for  all  the  Scandinavian  stations.  The  reduced  Q0  values 
resuited  because  the  Scandinavian  mb(Lg)  values  were  less  than  average. 

The  finite-difference  results  in  Figures  39  and  40  also  show  negatively  correlated  P  and 
SV  energy.  This  provides  a  preliminary  explanation  of  the  success  of  the  unified  yield  estima¬ 
tor  (U.S.  Congress,  Oflice  of  Technology  Assessment,  i 9 8 8 ) .  It  seems  that  measuring  all  pos¬ 
sible  phases  reduces  the  effects  of  uneven  energy  release  on  source  size  estimation.  To  under¬ 
stand  this  issue  in  a  more  quantitative  manner,  and  to  derive  an  optimal  weighting  scheme  to 
combine  all  phases,  theoretical  studies  with  numerical  simulations  are  necessary. 
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FIGURE.  CAPTIONS 

Figure  1.  Typical  geological  models  used  in  the  finite-difference  simulations  in  this  report, 
which  included:  (A)  model  5:  a  2  km  thick  self-similar  layer  of  mean  a  =  5  km/sec  and  10% 
variation,  superimposed  on  a  homogeneous  half  space  with  a  -=  6  km/sec.  (B)  A  2  km  thick 
Gaussian  layer  of  mean  ot  =  5  km/sec  and  10%  variation,  mean  grain  size  -  1  km.  (C)  model 
7:  a  3  km  thick  sedimentary  rock  of  toided  layers  with  sinusoidal  shape.  The  folded  layers 
have  wavelength  2  km,  peak-to-peak  amplitude  2.5  km.  Bach  iayer  is  uniform  in  material  pro¬ 
perty,  and  die  profile  of  a  is  Gaussian  distributed  with  mean  5  km/sec,  v  =  10%.  The  interface 
between  the  half  space  and  the  folded  layers  is  smoothed  (D)  model  9:  similar  to  (C)  except 
that  the  folded  layers  have  wavelength  5  km. 

Figure  2.  (E)  model  3:  a  rough  topography  suixrimposed  on  a  2  km  thick  Gaussian  layer. 
(F)  Same  as  (E)  except  that  the  layer  is  uniform  with  a  -  5  km/sec.  (G)  model  2:  similar  to 
(E),  except  for  self-similar  velocity  variation  and  different  topography.  (H)  model  4:  same  as 
(G;  except  for  uniform  layer.  The  rough  topography  was  generated  by  a  Markov  chain  with 
larger  transition  probability,  which  yielded  'more  frequent''  small-scale  elevation  changes. 

Figure  3.  Model  2.  P  wave  in  a  half  space  io  6.0  km/s,  p  =  3.55  km/s)  incident  upon  a  2 
km  layer  with  average  P-wave  velocity  of  5  km/s  and  a  self-similar  10%  rms  velocity  variation 
superimposed  by  a  gentle  topography  (indicated  in  the  0  sec  frame).  The  S-wave  velocity  is 
assumed  to  be  proportional  to  the  P-wave  velocity.  Darkness  of  the  snapshots  are  proportional 
to  the  displacement  amplitude.  Snapshots  of  the  displacement  field  are  shown  at  I  second 
intervals.  The  dilatational  strain  is  recorded  at  32  locations  at  a  depth  of  0.5  km  in  order  to 
infer  the  excitation  of  far-field  P  waves  from  explosion  sources.  Although  absorbing  boundary 
conditions  are  used,  care  must  be  taken  to  avoid  residual  reflections  from  the  sides  of  the  grid. 

Figure  4.  Model  7.  Same  as  Figure  3  except  that  a  3  km  thick  sedimentary  rock  of  folded 
layers  with  sinusoidal  shape.  The  folded  layers  have  wavelength  2  km,  peak-to-peak  amplitude 
2.5  kru.  Each  layer  is  uniform  in  material  property,  and  the  profile  of  a  is  Gaussian  distributed 
with  mean  5  km/sec,  v  =  10%.  The  interface  between  the  half  space  and  the  folded  layers  is 
smoothed  (see  Figure  1  (C)). 

Figure  5.  Model  8.  Same  as  Figure  3  except  the  layer  has  P-wave  velocity  of  5  km/s  and  a 
self-similar  20%  rms  velocity  variation  (indicated  in  the  0  sec  frame). 

Figure  6.  Model  14.  Same  as  Figure  5  except  that  the  gently  dipping  layers  totaling  3  km 
thick  are  superimposed  on  the  half  space.  The  sedimentary  layers  have  Gaussian  distributed 
profile  with  mean  5  km/sec,  v  =  10%. 

Figure  7.  Model  2.  S  wave  in  a  half  space  tot  =  6.0  km/s,  (3  =  3.55  km/s)  incident  at  52° 
upon  a  2  km  self-similar  layer  with  average  P-wave  velocity  of  5  km/s  and  10%  rms  velocity 
variation  superimposed  bv  a  gentle  topography  (indicated  in  the  0  sec  frame).  The  dilatational 
strain  is  recorded  at  32  locations  at  a  depth  of  0.5  km  in  order  to  infer  the  excitation  of  far- 
field  S  waves  from  explosion  sources. 

Figure  8.  Model  5.  Same  as  Figure  7  except  that  the  self-similar  layer  is  flatfindicated  in  the 
0  sec  frame). 
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Figure  9.  Model  9.  Same  as  Figure  7  except  that  the  3  km  thick  low  velocity  layer  consists  of 
lightly  folded  (5  km  wavelength,  2.5  km  amplitude)  0.5  km  thick  layers  with  10%  rms  velocity 
variation  with  respect  to  an  average  velocity  of  5  km/s. 

Figure  10  Model  14.  Same  as  Figure  7  except  that  the  half  space  is  superimposed  by  gently 
dipping  layers  totaling  3  km  thick.  The  sedimentary  layers  have  Gaussian  distributed  profile 
with  mean  5  km/sec,  v  =  10%. 

Figure  11.  Model  1.  Synthetic  far-iield  P-  (top)  and  SV-wave  (bottom)  inferred  by  the  princi¬ 
ple  of  reciprocity.  The  original  dilatational  strain  history  (5  Hz  low-pass)  responding  to 
incident  broadband  P  or  SV  plane  wave  recorded  at  32  locations  at  0.5  km  depth  in  the  refer¬ 
ence  model.  Shown  here  are  the  deconvolved  synthetics  corresponding  to  VSB  50  KT  in  hard 
rock.  The  peak  amplitude  of  these  synthetics  was  measured  and  compared  to  the  average  peak 
amplitude  of  the  reference  model. 

Figure  12.  Same  as  Figure  1 1  except  for  model  2. 

Figure  13.  Same  as  Figure  1 1  except  for  model  3. 

Figure  14.  Same  as  Figure  1 1  except  for  model  4. 

Figure  15.  Same  as  Figure  11  except  for  model  5. 

Figure  16.  Same  as  Figure  1 1  except  for  model  6. 

Figure  17.  Same  as  Figure  1 1  except  for  model  7. 

Figure  18.  Same  as  Figure  11  except  for  model  8. 

Figure  19.  Same  as  Figure  11  except  for  model  9. 

Figure  20.  Same  as  Figure  1 1  except  for  model  10. 

Figure  21.  Same  as  Figure  1 1  except  for  model  11. 

Figure  22.  Same  as  Figure  1 1  except  for  model  12. 

Figure  23.  San  e  as  Figure  11  except  for  model  13. 


Figure  24.  Same  as  Figure  1 1  except  for  model  14. 


P  [  Lgt 

figure  25.  Average  spectral  ratio  as  a  function  of  frequency.  Log{—~)  (upper)  and  Log{~ — ) 

•o  Lf>o 

(lower),  of  the  Model  1  relative  to  the  reference  model.  P  wave  response  of  model  1  in  the  0.5 
to  1.0  Hz  range  is  deficient  with  respect  to  the  reference  model  by  0.348  log  units,  while  the  S 
wave  response  is  0.063  log  unit  above  the  reference  model.  Vertical  bars  represent  the  stan¬ 
dard  error  of  a  single  observation. 


Figure  26.  Same  as  Figure  25  except  for  mode!  2. 
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Figure  27.  Same  as  Figure  25  except  for  model  3. 

Figure  28.  Same  as  Figure  25  except  for  model  4. 

Figure  29.  Same  as  Figure  25  except  for  model  5. 

Figure  30.  Same  as  Figure  25  except  for  model  6. 

Figure  31.  Same  as  Figure  25  except  for  model  7. 

Figure  32.  Same  as  Figure  25  except  for  model  8. 

Figure  33.  Same  as  Figure  25  except  for  model  9. 

Figure  34.  Same  as  Figure  25  except  for  model  10. 

Figure  35.  Same  as  Figure  25  except  for  model  1 1. 

Figure  36.  Same  as  Figure  25  except  for  model  12. 

Figure  37.  Same  as  Figure  25  except  for  model  13. 

Figure  38.  Same  as  Figure  25  except  for  model  14. 

Figure  39.  Averaged  P  and  SV(Lg)  peak  amplitudes  of  array  of  shallow  explosions  in  four¬ 
teen  crustal  models  sorted  with  the  log10(Lg/Lg0)  values.  Several  observations  are  obvious:  (1) 
Dipping  layers  (models  1 1  through  14)  generate  smaller  I.g  than  the  normalizing  model,  while 
they  all  generate  more  P  than  the  reference  model.  (2)  Media  with  topography  (e.g.  models  1 
through  4)  which  represent  CF.KTS  all  generate  more  Lg  than  the  normalizing  model,  while 
they  excite  less  P  due  to  strong  P  to  S  conversion.  (3)  Dipping  layers  (models  11  through  14) 
are  more  efficient  than  all  other  models  for  P  excitation.  Thus  mb{P)  and  mb{Lg)  appear  to  be 
negatively  correlated. 

Figure  40.  Same  as  Figure  39  except  the  P  and  SV(Lg)  excitations  are  measured  with  the 
averaged  spectral  content  in  [0.5,1]  Hz  band.  Crustal  models  with  topography  generate  more 
Lg  and  less  P  than  models  with  dipping  layers,  same  as  the  result  derived  from  peak  amplitude 
measurement. 

Figure  41.  Short  period  seismograms  of  two  Shagan  events  87I7I  (78.74E,  49.91N,  mb=6.1) 
and  87347  (78.85E,  49.96N,  mb=6.1),  and  a  Degelen  event  87198  (78.1  IE,  49.80N,  mb=5.8) 
recorded  at  CDSN  station  WMQ.  Each  trace  is  scaled  by  the  peak  amplitude.  Note  the  rela¬ 
tively  less  P  energy  (with  respect  to  Lg  energy)  in  the  Degelen  event  87198  as  compared  to 
Shagan  events  of  similar  magnitudes. 
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University  of  California 
Berkeley,  CA  9*720 

Professor  Thomas  H  .  Jordan 
Department  of  Earth,  Atmospheric 
and  Planetary  Sciences 
Mass  Institute  of  Technology 
Cambridge,  MA  02139 

Dr  .  Alan  Kafka 
Department  of  Geology  & 

Geophysics 
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Institute  of  Geophysics 
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Geosciences  Department 
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Department  of  Geosciences 
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University  of  Arizona 
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Cbrncll  University 
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Dr.  N  •  Biswas 
Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

Dr  .  G  .  A  .  Bollinger 
Department  of  Geological  Sciences 
Virginia  Polytechnical  Institute 
21044  Derring  Hall 
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Dr  .  James  Bulau 

Rockwell  Int’l  Science  Chnter 

1049  Cimino  Dos  Rios 

p  .0  .  Box  1085 

Thousand  Oaks,  CA  91360 

Mr  .  Roy  Burger 
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Schenectady,  N't  12309 

Dr.  Robert  Burrldgc 
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Professor  Jon  F.  Ciaorbout 
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Dr  •  Donald  Forsyth 
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Dr.  Anthony  Gnngi 
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Department  of  Geophysics 
Cbllcge  Station,  TX  77843 

Dr  •  Freeman  Gilbert 
Institute  of  Geophysics  & 

Planetary  Physics 
Univ .  of  Cilifornia,  San  Diego 
P  .0  .  Box  109 
La  Jolla,  CA  92037 

Mr  .  Edward  Gillcr 

Pacific  Seirra  Research  Cbrp  . 

1401  Wilson  Boulevard 
Arlington,  VA  22209 

Dr.  Jeffrey  W.  Given 
Sierra  Geophysics 
11255  Kirkland  Way 
Kirkland,  WA  98033 

Dr  .  Arthur  Lerncr-Lam 

I.nmont-Dohor ty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 

Dr  .  L  .  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  CA  30332 
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La  Jolla,  CA  92093 
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SRI  International 
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Cbrncll  University 
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Dr.  Alan  S.  Ryall,  Jr. 

G'nter  of  Seismic  Studies 
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Dr  .  David  G  .  Simpson 
Lamont-Dohcrty  Geological  Observ  . 

of  Gblumbia  University 
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University  of  Utah 
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Geophysics  Program 
University  of  Washington 
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Dr.  Jerry  Chrter,  Dr.  Paul  Pomeroy 
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Dr  .  L  •  Sykes 

Lamont  Doherty  Geological  Observ  . 
Cblumbia  University 
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Dr  .  Pradoep  Taiwan! 

Department  of  Geological  Sciences 
University  of  South  Carolina 
Cblumbia,  SC  29208 

Dr  .  R  •  B  .  Tittmann 

Rockwell  International  Science  Ctenter 
1049  Chmino  Dos  Rios 
P  .0  .  Box  1085 
Thousand  Oaks,  CA  91360 

Weidlinger  Associates 
ATTN:  Dr.  Gregory  Wojcik 
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Harvard  University 
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Dr .  Fetor  Basham 
Earth  Physics  Branch 
Geological  Survey  of  Canada 
1  Observatory  Crescent 
Ottawa,  Ontario 
CANADA  K1A  OY3 

Dr  .  Eduard  Berg 
Institute  of  Geophysics 
University  of  Hawaii 
Honolulu ,  III  96822 

Dr.  Michel  Bouchon  -  Universite 

Sc  lent l f Ique  et  Medlcalc  de  Grenob 
Lab  de  Geophysique  -  Interne  et 
Tectonophysique  -  1  .R  .1  .G  .M-B  .P  . 
38402  St.  Martin  D'lleres 
Ccdcx  IRAN  CE 

Dr.  Hilmar  Bungum/NTN F/NORS AR 
P -0 .  Box  51 

Norwegian  Cbuncil  of  Science, 

Industry  and  Research,  NORSAR 
N-2007  Kjeller,  NORWAY 

Dr .  Michel  Chmpillo 
I  .R.I  .C  .M.-B.P.  68 
38402  St  .  Martin  D'Heres 
Chdex,  IRANCE 

Dr.  Kin-Y.ip  Chun 
Geophysics  Division 
Physic.'.'  Department 
University  of  Toronto 
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Dr.  Alan  Douglas 
Ministry  of  Defense 
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UNITED  KINGDOM 
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Fed.  Inst.  For  Geosciences  A  Nat  *1  Res. 
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FEDERAL  REPUBLIC  OF  GERMANY 
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NTNF/  NORSAR 
r  .0  .  Box  51 

N-2007  Kjeller,  NORWAY 
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Mr.  Peter  Marshall,  Procurement 
Executive,  Ministry  of  Defense 
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UNITED  KINGDOM  (3  copies) 

Dr  .  Ben  Menahelm 

Weizman  Institute  of  Science 

Rehovot,  ISRAEL  951729 

Dr.  Svcin  Mykkeltveit 
NTN F/NORSAR 
P  .0  .  Box  51 

N-2007  Kjeller,  NORWAY  (3  copies) 

Dr  .  Robert  North 
Geophysics  Division 
Geological  Survey  of  Chnada 
1  Observatory  crescent 
Ottawa,  Ontario 
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Dr  .  Frode  Ringdal 
NTN F/NORSAR 
P  .0  .  Box  51 

N-2007  Kjeller,  NORWAY 
Dr.  Jorg  Schlittcnhardt 

Federal  Inst,  for  Geosciences  &  Nat’l  Res. 

Postfach  510153 

D-3000  Hannover  51 

FEDERAL  REPUBLIC  OF  GERMANY 

University  of  Hawaii 
Institute  of  Geophysics 
ATTN:  Dr.  Daniel  Walker 
Honolulu,  HI  96822 
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Dr.  Ramon  Cibrc,  S  .J  . 
c/o  Hr .  Ralph  Buck 
Economic  (bnsular 
American  Embassy 
ATO  Miami,  Florida  34032 

Professor  Peter  Harjes 
Institute  for  Geophysik 
Rhur  Univcrsity/Bochum 
P  .0 .  Box  102148,  4630  Bochum  I 
FEDERAL  REPUBLIC  OF  GERMANY 

Professor  Brian  L  .N  •  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
G  .P  .0  ■  Box  4 
Chnbcrra  2601 
AUSTRALIA 

Dr  .  B  .  Hass i non 

Sociote  Radiomana 

27,  Rue  Claude  Bernard 

7,005,  raris,  FRANCE  (2  copies) 

Dr .  Pierre  Mochler 
Societe  Radiomana 
27 ,  Rue  Qaude  Bernard 
75005,  Paris,  FRANCE 
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Dr  .  Ralph  Alewine  III 

darpa/nmro 

1400  Wilson  Boulevard 
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Dr  •  Peter  Basham 
Geological  Survey  of  Canada 
1  Observatory  Creseut 
Ottovn,  Ontario 
CANADA  K1A  0Y3 

Dr  •  Robert  Blandford 
DARPA/NMRO 

1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Sandia  National  Laboratory 
ATTN:  Dr.  H.  B.  Durham 
Albuquerque,  NM  87185 

Dr  •  Jack  Evcrnden 
USGS -Earthquake  Studies 
345  Middlcf ield  Road 
Henlo  Park,  CA  94025 

U  .S  .  Geological  Survey 
ATTN  :  Dr  .  T  .  Hanks 
Nat'l  Earthquake  Resch  Chnter 
345  Middlef ield  Road 
Menlo  Park,  CA  94025 

Dr  .  James  Hannon 

Lawrence  Livermore  Nat'l  Lab. 

P  .0  .  Box  808 
Livermore,  CA  94550 

U  .S  .  Arms  Obntrol  &  Disarm.  Agency 

ATTN:  Mrs.  M.  lloinkcs 

Div  .  of  Multilateral  Affairs 

Room  5499 

Washington,  D-C.  20451 

Paul  Johnson 

ESS-4 ,  Mill  Stop  J979 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87545 
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&  Engineering 
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Menlo  Park,  CA  94025 

Dr  .  William  Leith 

uses 

Mail  Stop  928 
Reston,  VA  22092 

Dr.  Robert  Masse' 

Box  25046,  Mail  Stop  967 
Denver  Federal  Genter 
Denver,  Cblorado  80225 

Dr.  Keith  K.  Nakanishi 

Lawrence  Livermore  National  Laboratory 

F  .0  .  Box  008,  L-205 

Livermore,  CA  94550  (2  copies) 

Dr  .  Carl  Newton 

Los  Alamos  National  Lab . 

I’  .0  .  Box  1663 

Mall  Stop  035,  Group  E553 
Los  Alamos,  NM  87545 

Dr.  Kenneth  H.  Olsen 

Los  Alamos  Scientific  Lab  • 

Pont  Office  Box  1663 
Los  Alamos,  NM  87545 

Howard  J  .  Patton 

Lawrence  Livermore  National  Laboratory 
P.0.  Box  808,  L-205 
Livermore,  CA  94550 

HQ  AFTAC/T7- 

Attn:  Dr.  Frank  F.  Pilotte 
Patrick  AFB,  Florida  32925-6001 

Mr  .  Jack  Rnchlin 
USGS  -  Geology,  Rm  3  C136 
Mall  Stop  928  National  Cbnter 
Reston,  VA  22092 

Robert  Rei nke 
AFWL/NTESG 

Kirtland  AFB,  NM  87117-6008 
HQ  AFTAC/TGR 

Attn:  Dr.  George  II .  Rothe 
Patrick  AFB,  Florida  32925-6001 
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Donald  L.  Springer 

Lawrence  Livermore  National  Laboratory 
P.0.  Box  808,  L-205 
Livermore,  CA  94550 

Dr  •  Lawrence  Turnbull 
OSWR/NED 

Central  Intelligence  Agency 
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Washington,  D.C.  20505 

Dr  .  Thomas  Weaver 
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Office  of  the  Secretary  Defense 

DDR  4  E 
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ATTN:  Technical  Library 

Washington,  DC  20305 

Director,  Technical  Information 
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Defense  Technical 
Information  Center 
Cimeron  Station 
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Washington,  D.C.  20301 
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